Stretched exponential relaxation for growing interfaces in quenched disordered media 
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We study the relaxation for growing interfaces in quenched disordered media. We use a directed 
percolation depinning model introduced by Tang and Leschhorn for 1 + 1-dimensions. We define the 
two-time autocorrelation function of the interface height C(t',t) and its Fourier transform. These 
functions depend on the difference of times t—t' for long enough times, this is the steady-state regime. 
We find a two-step relaxation decay in this regime. The long time tail can be fitted by a stretched 
exponential relaxation function. The relaxation time r a is proportional to the characteristic distance 
of the clusters of pinning cells in the direction parallel to the interface and it diverges as a power law. 
The two-step relaxation is lost at a given wave length of the Fourier transform, which is proportional 
to the characteristic distance of the clusters of pinning cells in the direction perpendicular to the 
interface. The stretched exponential relaxation is caused by the existence of clusters of pinning cells 
and it is a direct consequence of the quenched noise. 
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I. INTRODUCTION 

For decades the investigation of growing surfaces and 
interfaces has attracted much attention due to its impor- 
tance in many fields, such as motion of liquids in porous 
media, growth of bacterial colonies, crystal growth, fronts 
of fire, etc. In these problems we have a nonequilib- 
rium interface. The d-dimensional interface described 
by a single valued- function h(x, t) evolves in a d + 1- 
dimensional medium. The disorder affects the motion of 
the interface and leads to its roughness. A phenomeno- 
logical non-linear Langevin equation, the Kardar-Parisi- 
Zhang equation (KPZ) [1], and the directed percolation 
depinning (DPD) models [2,3] have been used in order 
to study growing interfaces. Two main kinds of disor- 
der have been proposed in these models: the annealed 
noise that depends only on time and the quenched disor- 
der due to the inhomogeneity of the media, which does 
not depend on time. In the DPD models the disorder is 
quenched and they describe very well some experiments 
such as the growth of bacterial colonies and the motion 
of liquids in porous media. These models were proposed 
simultaneously by Tang and Leschhorn [2] and Buldyrev 
et al [3]. 

In many glassy systems a non-exponential relaxation 
is found when they are close to some temperature above 
to the static transition. As an example, in structural 
glasses a two-step relaxation decay is found near the so 
called "ideal glass transition" [4]. The long relaxation 
step has the stretched exponential form 



f(t) = f e X p[-(t/ Ta f] 



(1) 



where < (3 < 1 does not depend on the tempera- 
ture. There are two mechanisms driving non-exponential 
relaxation. In disordered systems such as spin glasses 



that behavior is caused by the existence of non-frustrated 
ferromagnetic-type clusters of interactions [5] which is a 
direct consequence of the quenched disorder [6] . Another 
mechanism in frustrated systems is based on the perco- 
lation transition of the Kasteleyn-Fortuin and Coniglio- 
Klein cluster [7], here disorder is not needed to obtain 
non-exponential relaxation [8]. Recently, Colaiori and 
Moore [9] have found a stretched exponential relaxation 
for the KPZ equation with annealed noise. 

In this paper, we use the DPD model proposed by 
Tang and Leschhorn (TL) [2] in order to investigate the 
relaxation of the two-time autocorrelation functions in 
quenched disordered media at the steady-state regime. 
We relate the relaxation properties to the clusters of pin- 
ning cells. The paper is organized as follows. In Sec. II 
we present the model and some properties of the clusters 
of pinning cells. In Sec. Ill the steady-state relaxation is 
studied. Finally, in Sec IV we present some conclusions. 



II. THE MODEL 

In the TL model for 1 + 1-dimensions [2] , the advance 
of the fluid through the media is modeled by a driving 
force p, while the disorder of the media, that brakes this 
advance, is represented by a quenched noise in the sub- 
stratum. The interface grows in a square lattice of edge L 
with periodic boundary conditions. We assign a random 
pinning force g(r) uniformly distributed in the interval 
[0, 1] to every cell of the square lattice. For a given ap- 
plied driving force p > , we can divide the cells into 
two groups: those with g(r) < p (free cells), and those 
with g(r) > p (pinning cells). Denoting by q the den- 
sity of pinning cells on the lattice, we have q = 1 — p 
for < p < 1 and q = for p > 1. The interface is 
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specified completely by a set of integer column heights 
hi (i = 1, . . . , L). At t = all columns are assumed 
to have the same height, which is zero. During growth, 
a column is selected at random, say column i, and its 
height is compared with those of the neighbor columns 
(i — 1) and (i + 1). The growth event is defined as follows. 
If hi is greater than either or hi+i by two or more 
units, the height of the lower of the two columns (i — 1) 
and (i + 1) is incremented in one (in case of the two being 
equal, one of them is chosen with equal probability). In 
the opposite case, hi < min(/ij_i, h i+ i) + 2, the column 
i advances by one unit provided that the cell to be oc- 
cupied is a free cell. Otherwise no growth takes place. 
In this model, the time unit is defined as one growth at- 
tempt. In numerical simulations at each growth attempt 
the time t is increased by St, where St — 1/L. Thus, after 
L growth attempts the time is increased in one unit. In 
our simulations we use L = 10000 and take the averages 
over 100 different realizations of quenched noise. 

A. Clusters of pinning cells 

As it has been shown in Ref. [2], this model has a de- 
pinning transition at a driving force p c = 0.461. For 
driving forces below the critical one p c , the advance of 
the interface is halted (pinning phase), while above this 
driving force the interface moves without stopping (mov- 
ing phase). At the transition, the characteristic length 
£ of the pinned regions diverges. A directed percolation 
cluster of pinning cells which extends over the whole sys- 
tem appears in the pinning phase. In the moving phase, 
a typical connected cluster of pinned cells extends over a 
distance of the order of £|| in the direction parallel to the 
interface and a distance of the order of £j_ in the direc- 
tion perpendicular to the interface. On both sides of the 
percolation transition, the two lengths have a power-law 
behavior £|| ~ \p — p c \~" 11 and £j_ ~ \p — p c \^ v± , with 
u\\ = 1.733 ± 0.001 and v ± = 1.097 ± 0.001. £_l sets a 
characteristic scale for the height while £|| sets character- 
istic scales for both the distance parallel to the interface 
and the time. For the mean interface height the scaling 
form H(t) ~ £±®{t/{;\\) is obtained, denoting $ a scaling 
function which is different for the two phases. In the mov- 
ing phase there is a crossover from a power-law growth 
H{t) ~ t v± / v W at t < £|| to a linear behavior H(t) = vt 
at t » £|| . The steady-state velocity can be expressed as 
v(p) ~ (p-p c ) v ii 

III. STEADY-STATE RELAXATION 

We define the two-time autocorrelation function of the 
surface height as 

c(t',t) = ^s[h j (t')-h j (t)} , (2) 

3 



where S [x] is the delta function. Its Fourier transform is 

C l (i',() = iVe-[' , '( f '»-'''l t »l l , (3) 
Li * 

3 

where k is the wave number. 
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FIG. 1. C k (t - t') for p = 0.5 and k = ir and for different 
initial times t' = 10 2 (long dashed line), 10 3 (dashed line), 
10 4 (dotted line), and 10 5 (solid line). Dotted and solid lines 
overlap for any t — t'. 

For long enough times, t' 3> £y , these functions depend 
only on the difference of times t — t' , this is the steady- 
state regime where H (t) = vt. This regime is reached 
at longer times when we approach to the critical driving 
force p c . Fig. 1 shows Ck{t — t') for p — 0.5 and k = ir 
and for different initial times t' = 10 2 , 10 3 , 10 4 , and 10 5 . 
As we see Ck(t — t') is independent of t' when t' > 10 4 
for this value of p. 
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FIG. 2. Ck(t) in the steady-state regime for p — 0.5 and 
k = 7r (solid line), tt/2 (dotted line), it/ 3 (dashed line), 7r/4 
(long dashed line), and it/ 5 (dot-dashed line). 

We find a two-step relaxation decay in the steady-state 
regime. We can see in Fig. 2 that the time interval of 
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the first and second relaxation step depends on the wave 
number k. Nevertheless, the form of the second relax- 
ation step does not depend on k. For small enough k we 
only have one step relaxation process. So, there is a wave 
number k e where the two-step relaxation decay is lost for 
k < k e . Fig. 3 shows Ck{t) for k = tt and different values 
of the driving force p. The two-step relaxation decay is 
observed from the highest value of p — 0.95, but the time 
interval of the second step increases when p is decreased, 
i.e. when the system approaches to the criticality. This 
behavior is also found in other glassy systems, where the 
time interval of the second relaxation step increases when 
the systems approach to the critical temperature. As we 
can see in Fig. 4, the second relaxation step can be fit- 
ted by a stretched exponential relaxation function Eq. 
(1) with the exponent (3 — 0.805 ±0.05. This exponent is 
in practice independent of p and it brings to the master 
equation 

Ck(t) = C k (t/r a ) (4) 

for t > T a , where Ck(t/r a ) does not depend on p. In 
glassy systems Eq. (4) is also called time-temperature su- 
perposition principle [4], because the temperature plays 
the role of the driving force in that systems. In the in- 
set of Fig. 4 we show an equivalent time-driving force 
superposition principle for our system. This stretched 
exponential relaxation means that in the system there is 
a broad distribution of relaxation times [10]. 
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FIG. 3. Ck(t) in the steady-state regime for k — tt and 
p = 0.95, 0.9, 0.85, 0.8, 0.75, 0.7, 0.65, 0.6, 0.68, 0.56, 0.55, 
0.54, 0.53, 0.52, 0.51, 0.5, 0.49, 0.48, 0.475, and 0.47 (from 
left to right). 

The relaxation time r Q can be obtained from the fit of 
Ck(t) with a stretched exponential function, it is shown 
in Fig. 5. We see that it is very well fitted by a power 
law r a oc (p — Pc) _I/ " where p c = 0.462 ± 0.001 and vu = 
1.733 ± 0.001. This means that r a is proportional to 
the characteristic distance of the clusters of pinning cells 
in the direction parallel to the interface £y, which also 



diverges as r Q ~ £|| ~ (p — p c ) "n . 
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FIG. 4. Log- log plot of Ck(t) where k = tt and p = 0.56, 
0.55, 0.54, 0.53, 0.52, 0.51, 0.5, 0.49, 0.48, 0.475, and 0.47 
(from left to right). Dashed curves are fitting functions cor- 
responding to the stretched exponential functions. Inset: 
time-driving force superposition principle. The dashed curve 
is a stretched exponential function with f3 = 0.8. 

We can obtain a wave length A e = 7r/fe e in the direction 
perpendicular to the interface from the wave number k e 
where the two-step relaxation is lost. In the inset of 
Fig. 5 we show A e as a function of (p — p c ). This length 
diverges as a power law A e oc (p — p c )~ u± with p c — 
0.46±0.01 and v± = l.lzbO.Ol. So that, A e is proportional 
to the characteristic distance of the clusters of pinning 
cells in the direction perpendicular to the interface, A e ~ 
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FIG. 5. Log-log plot of the relaxation time r a , obtained 
by the stretched exponential fit in Fig. 4, as a function of 
p -p c , for p = 0.56, 0.55, 0.54, 0.53, 0.52, 0.51, 0.5, 0.49, 
0.48, 0.475, and 0.47. The solid curve is a power law function 
r a = 0.275(p - 0.462)- 1 ' 733 . Inset: log-log plot of A e as a 
function of p - p c , for p = 0.54, 0.53, 0.52, 0.51, 0.5, 0.49, 
0.48, 0.475, and 0.47. The solid curve is a power law function 
A e = 1.25(p -0.46) -1 - 1 . 
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We see that the characteristic length of the clusters of 
pinning cells in the direction parallel to the interface £y 
sets the time scale of the stretched relaxation function 
Cfc(V£||)- On the other hand, the stretched relaxation 
step is lost for A > A e , that is for A > So, the 

stretched exponential relaxation is caused by the clusters 
of pinning cells. 

IV. CONCLUSIONS 

We have studied relaxation properties for growing in- 
terfaces in quenched disordered media. We have used 
the TL model in which properties of clusters of pinning 
cells arc known. We have studied the relaxation prop- 
erties of the Fourier transform of the autocorrelation of 
the surface height and found a two-step relaxation pro- 
cess in which the second step is well fitted by a stretched 
relaxation function with (3 = 0.805 ±005. The relaxation 
time diverges as a power law and it is proportional to the 
characteristic distance of the clusters of pinning cells in 
the direction parallel to the interface. The form of the 
second step relaxation does not depend on the wave num- 
ber of the Fourier transform. This step is lost for a given 
wave length A e which is proportional to the characteristic 
distance of the clusters of pinning cells in the direction 
perpendicular to the interface. From these results, we 
can say that the stretched exponential relaxation behav- 
ior is caused by the clusters of pinning cells, which is a 
direct consequence of the quenched noise as it happens 
in other glassy systems [6] . 



[6] A. Ficrro, A. de Candia, and A. Coniglio, Phys. Rev. E 

56, 4990 (1997). 
[7] CM. Fortuin and P.W. Kasteleyn, Physica (Amaster- 

dam) 57, 536 (1972); A. Coniglio and W. Klein, J. Phys. 

A 12, 2775 (1980). 
[8] A. Fierro, G. Franzese, A. de Candia, and A. Coniglio, 

Phys. Rev. E 59, 60 (1999). 
[9] F. Colaiori and M.A. Moore, Phys. Rev. E 63, 057103 

(2001). 

[10] See, e.g., R. Richert and A. Blumen, in Disorder Effects 
on Relaxational Processes, edited by R. Richert and A. 
Blumen (Spinger-Verlag, Berlin, 1994), p. 1. 



ACKNOWLEDGMENTS 



This work was supported in part by the project No. 
PI-60/00858/FS/01 from the Fundacion Seneca, Region 
de Murcia. 



[1] M. Kardar , G. Parisi, and Y.C. Zhang, Phys. Rev. Lett. 

56, 889 (1986). 
[2] L.H. Tang and H. Leschhorn, Phys. Rev. A 45, R8309 

(1992). 

[3] S.V. Buldyrcv, A.L. Barabasi, F. Caserta, S. Havlin, H.E. 
Stanley, and T. Vicsek, Phys. Rev. A 45, R8313 (1992). 

[4] W. Gotze, in Liquids, Freezing and Glass Transition, 
edited by J. P. Hansen, D. Levesque, and P. Zinn- Justin 
(Elsevier, Amsterdam, 1991). 

[5] M. Randeria, J. P. Sethna, and R.G. Palmer, Phys. Rev. 
Lett. 54, 1321 (1985); F. Cesi, C. Maes, and f. Martinelli, 
Commun. Math. Phys. 188, 135 (1997); G. Franzese and 
A. Coniglio, Phys. Rev. E 58, 2753 (1998). 



4 



